A brown dwarf candidate in the Praesepe Open Cluster 
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ABSTRACT 

^ ■ We present optical and infrared observations of RPrl, a faint (/ = 21.01) and very red 

OO . object {I — K = 4.57) discovered in a deep CCD survey, covering an area of 800 square arcmin 

of the Praesepe open cluster. A low resolution spectrum shows that RPrl is a very late object, 
■ the latest object in Praesepe for which a spectrum has been taken to date. Our estimates give a 

, mass between 0.063 and 0.084 A/q, and indicate that RPrl may turn out to be the first brown 

' dwarf in this cluster. 

Subject headings: open clusters and associations: individual (Praesepe) — stars: low-mass, 
O ' brown dwarfs 



Introduction 



Young, close, rich clusters make ideal hunting grounds for brown dwarfs (BDs). The advantage of 
' searching in these clusters is that important parameters like age, distance, extinction, metallicity are 

relatively well constrained. Moreover, the probability of detection will be increased by the higher intrinsic 
luminosity of the cluster substellar objects and the higher stellar density in the cluster with respect to the 
field. To date, searches for substellar bodies have been carried out quite extensively in the Pleiades cluster 
(see Hambly 1998 for a review), where recently BDs down to 0.045Mq have been confirmed (Zapatero 
Osorio et al. 1997a) and many candidates down to O.OBOMq wait for confirmation. 

A comparison of the lower end of the mass function of the Pleiades with that in other open clusters 
would be very valuable in order to track the evolution of very low mass objects around the substellar limit 
at different ages and set constraints on predictions by models. Therefore, an extension of the search for 
BDs to clusters of different ages will add to our understanding of evolution in the substellar regime. An 
interesting cluster besides the Pleiades is the Praesepe open cluster (M44). Several searches have already 
been performed on this cluster. The very low main sequence of Praesepe was traced by the photometric 
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and proper motion survey by Hambly et al. (1995), and its mass function was determined down to O.OSMq 
by Williams, Rieke & Stauffer (1995). Recently, Pinfield et al. (1997) published new BD candidates, with 
magnitudes down to / = 21.5. We discuss here low resolution spectroscopy and IR and optical photometry 
of an object that may turn out to be the first bona fide brown dwarf in Praesepe. 

2. Observations and data reduction 

The object presented in this paper was selected as a substellar candidate in a deep R, I CCD survey 
in the Praesepe open cluster conducted in February 1996 at the prime focus of the 2.5 m Isaac Newton 
Telescope, La Palma. The survey covered an area of 800 square arcmin, complete down to / = 21.2, 
R = 22.2. The survey sample, data reduction and photometric analysis will be described in a forthcoming 
paper. Instrumental magnitudes for RPrl were derived from aperture photometry, and conversion to the 
standard R, I Cousins system was performed observing standard stars from Landolt (1992). In agreement 
with the notation in Stauffer et al. (1989) and Rebolo, Zapatero Osorio, & Martin (1995), we named 
our object Roque Praesepe 1, hereafter referred to as RPrl, for brevity. The coordinates of RPrl are 
a = 8'^38™02';9, S = +19°25'46'.'7 (Equinox 1950, Epoch 1997.11). In Figure |l] we show a finding chart for 
RPr 1 (/ image) . 

Follow up infrared photometry was performed with the WHIRCAM camera at the 4.2 m William 
Herschel Telescope on October 30, 1996. The camera was equipped with a InSb 256 x 256 detector array, 
using the J and K' filters. The total integration time per filter was 200 s, each final image being the 
coaddition of 5 dithered exposures of 40 s each. Data were processed using standard techniques within 
the IRAfQ environment. Instrumental aperture magnitudes were corrected for atmospheric extinction and 
transformed into the UKIRT system by observing the standard star F16 (Casali & Hawarden 1992). The 
K' data were translated into K magnitudes using the conversion in Wainscoat & Cowie (1992). Such 
conversion makes use of the H — K color. We have no H magnitude for RPrl, so we adopted the average 
H — K between Calar 3 and Teide 1, the reddest objects in the Zapatero Osorio, Martin, & Rebolo (1997b) 
sample, obtaining H — K — 0.5 ± 0.2. In the Wainscoat & Cowie conversion our ±0.2 error on H — K 
translates into a ±0.04 uncertainty, which is included in the error on K quoted in Table ^, where the results 
of the optical and infrared photometry of RPr 1 are summarized. 

Table 1: Photometry and spectral type of RPr 1 



R 


/ 


J 


K 


Sp. Type 


23.5±0.2 


21.01±0.05 


17.7±0.1 


16.44±0.08 


- A/8.5 



Low-resolution optical spectroscopy with the red arm of the ISIS spectrograph was taken at the William 
Herschel Telescope on 1996 December 9. The detector was a TEK 1024 x 1024 CCD and the grating R158R 
gave a nominal dispersion of 2.9 A per pixel in the range 650-915 nm. Two spectra were taken (exposure 
times 3600 and 5300 s). Spectra were reduced by a standard IRAF procedure. The instrumental signature 
was removed making use of the standard G191B2B. 



^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under contract with the National Science Foundation. 
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3. Analysis 

In Figure ^, RPr 1 is located in the K ys. I ~ K color-magnitude diagram, together with infrared 
photometry from Hodgkin et al. (1998) for proper motion members fr-om Hambly et al. (1995) and objects 
from the survey of Pinfield at al. (1997). The solid lines are two isochrones from solar metallicity models by 
Baraffe et al. (1998), for 400 and 900 Myr, shifted to the distance of Praesepe, which we set equal to 177 pc 
(m — M — 6.24), the value determined by Hipparcos (Mermilliod et al. 1997). Although various ages have 
been reported for Praesepe (cf. Hambly et al. 1995), it is generally believed (see e.g. Stauffer 1996) that this 
cluster has similar age as the Hyades 600 Myr). A conservative estimate lies between 400 Myr, the value 
reported by Allen (1973), and 900 Myr, as reported by Vandenberg & Bridges (1984). Following Hambly 
et al. (1995) we will assume extinction and reddening to be zero in front of Praesepe. We note that RPr 1 
fits the photometric sequence defined by Hodgkin et al. (1998) extremely well. This is very promising as far 
as the confirmation of cluster membership is concerned. 

In Figure ^ we present the best of our two low-resolution optical spectra of RPr 1 . This spectrum 
(exposure time 5300 s) shows prominent TiO and VO molecular absorption bands, typical of very late M 
stars, and rather strong atomic lines of Na I (818.3-819.5 nm) and K I (766.5-769.9 nm), indicating a dwarf 
nature. In order to derive an accurate spectral type, we have measured the pseudocontinuum spectral ratios 
defined by Martin, Rebolo, & Zapatero Osorio (1996) and the VO spectral index, defined by Kirkpatrick, 
Henry, & Simons (1995). Of the five PC indices by Martin et al. (1996), PC5 is outside our spectral range, 
and PCI, PC2 lie in the noisiest region of our spectrum. Thus, only PC3 and PC4 have been considered. 
According to Martin et al.'s results these two indices are the most sensitive to the spectral type. We obtain 
PC3 = 1.80, PC4 = 3.06, VO = 1.13, which give a spectral type ranging between M8.2 and M9. 

The above spread may come from the poor signal-to-noise ratio of our spectrum. However, our 
measurements confirm that we are dealing with a very late object, the latest object in Praesepe for which a 
spectrum has been taken to date. In addition, the I — K oi RPr 1 (4.57) corresponds to those of M9 dwarfs, 
both in the field (Kirpatrick & McCarthy 1994) and in the Pleiades (Zapatero Osorio et al. 1997c). We 
have to remark that our derivation of spectral type could be affected by the Praesepe metallicity, slightly 
higher than solar (Fe/H = 0.14, Reglero & Fabregat 1991). Moreover, Jones & Tsuji (1997) argue that 
molecular indices should not be used for spectral classification of very late-type objects, due to possible 
grain condensation in the outer layers of their atmospheres. 

We do not detect any Ha emission in the spectrum of RPrl. However, the signal-to-noise ratio in the 
Ha region is very low, so we estimate an upper limit of about 3 A (Icr) for Ha equivalent width. Williams 
et al. (1994) measured Ha emission in low mass stars in Praesepe; for M3-M4 stars (the latest type in 
their sample) they found IVA(Ha) between 2.5 and 5 A. To our knowledge, no data are available for later 
types in this cluster. Later spectral types were investigated by Stauffer et al. (1997) in the Hyades; again, 
equivalent widths not exceeding 5 A were found. Our upper limit is in agreement with these results, and is 
consistent — or, at least, not in contrast — with RPr 1 membership to Praesepe. We note that in Roque 4, 
a Pleiades member with similar spectral type. Ha emission has not been observed (Zapatero Osorio et al. 
1997a). These authors point out that a membership criterion based on Ha emission might not be valid for 
spectral types as late as M9. 

An attempt to measure the radial velocity in the spectrum of Figure ^ by cross-correlation with 
templates observed using the same instrumental configuration gave a value of about zero with a la errorbar 
of ±40 km s~^. The radial velocity of the Praesepe cluster (-f34.3 km s"^, Mermilliod et al. 1990) lies in 
the large range defined by the above uncertainty, but obviously higher resolution spectroscopy and better 
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signal-to-noise ratio are required in order to confirm RPr 1 membership through radial velocity. 

From a comparison of the K apparent magnitude of RPrl with the Baraffe et al.'s (1998) isochrones, 
we can estimate the mass of this object. Considering a possible range of distances between m — M = 6.24 
(Mermillod et al. 1997) and 6.05 (Reglero & Fabregat 1991) and ages between 400 and 900 Myr, we get a 
mass range between 0.063 and 0.084 Mq. Note that the main uncertainty on the mass of RPr 1 comes from 
the wide range of possible ages of the cluster. 

An average of different temperature calibrations (Bessell & Stringfellow 1993; Kirkpatrick 1995; 
Jones et al. 1994; Brett 1995; Jones et al. 1996) gives an effective temperature of 2400^22° K for RPrl, 
considering that its spectral type is between M8.2 and M9. The luminosity of RPr 1 was determined by 
taking an average of the values obtained from three different calibrations: bolometric magnitude vs. I — K 
(Tinney et al. 1993), K bolometric correction vs. I — K (Tinney et al. 1993), and bolometric magnitude vs. 
effective temperature (Jones et al. 1994). The result is log{L/LQ) = —3.60 ±0.16. 



4. Discussion and conclusions 

Every point of the above analysis supports the membership of RPr 1 in the Praesepe cluster. However, 
we still cannot rule out that we might be dealing with a very late M dwarf of the field. A conclusive 
answer would come from the radial velocity. In absence of such measurement, further support to the 
RPr 1 membership can be provided by the knowledge of the probability for a late M field dwarf to appear 
disguised as RPr 1. Below we try to estimate how many M8.2-M9 field dwarfs with a magnitude equal to 
that of RPr 1 within the photometry errors {2a, to be conservative) can be expected in the area we surveyed 
(800 square arcmin). 

From the literature (Leggett 1992; Kirkpatrick & McCarthy 1994; Tinney et al. 1995; Monet et al. 
1992) we derive the absolute K magnitudes of M8 and M9 stars as 10.02 ±0.03 and 10.41 ±0.13, respectively. 
Thus, we assume that any contaminating M8.2-M9 dwarf would have 10.0 ^ Mk <, 10.4. Considering a 2a 
uncertainty around K = 16.44 we find that such field dwarf would lie between 150 and 208 pc. From results 
in Kirkpatrick et al. (1994) we estimate that the density of M8.5 field dwarfs in the solar neighborhood is 
0.0029 pc"'^. Taking into account a correction for the decrease in density when increasing the height above 
the galactic plane, with a scale height of 350 pc (Mihalas & Binney 1981), we obtain that the number of 
contaminating very late M field dwarfs that we expect to find in our survey with the same K magnitude as 
RPr 1 is about 0.24. In other words, we would have needed to survey an area about 4 times larger in order 
to find one field dwarf with the RPr 1 characteristics. 

On the other hand, we have to calculate how many BDs like RPr 1 we expect in the cluster. In the 
hypothesis that the mass function of Praesepe in the substellar regime follows the trend found by Williams 
et al. (1995) for masses down to 0.2 Mq (a power law with slope —1.34), we would expect 3 objects in the 
mass range 0.063-0.084 Mq in 800 square arcmin. Actually, in this calculation we have to use the mass 
range obtained adopting a unique age. If we use 600 Myr we obtain a mass range 0.074-0.080 Mq and 1 
object in 800 square arcmin. Even in these more restrictive conditions, RPr 1 has a probability 4 times 
higher to belong to Praesepe than to be a field star. 

The above calculations strongly support RPrl's membership. This object results the faintest and 
coolest member of Praesepe and may turn out to be the first bona fide brown dwarf in this cluster. The 
possible rise in the mass function of the Pleiades beyond the substellar frontier suggested by the data of 
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Zapatero Osorio ct al. (1997a, b) and our finding of RPr 1 encourage a deeper extension of the survey in a 
larger area of Praesepe. 
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Fig. 1. — Finding chart of RPr 1 (7 filter). Size is 2 x 2 arcmin. North is up and East at right. 

Fig. 2. — The K vs. I — K diagram for very low mass stars in Praesepe. Together with RPr 1 (big dot with 

errorbars), wc show two isochrones, at 400 (upper curve) and 900 Myr (lower curve), from Baraffe et al. 
(1998); data from Hodgkin et al. (1998) for proper motion members from Hambly et al. (1995; filled circles) 
and objects from Pinfield et al. (1997; open squares). In the isochrones, masses of 0.08 Mq (stars) and 0.07 
Mq (crosses) are indicated. Errorbars are at la. 

Fig. 3. — Low-resolution spectrum of RPr 1. Relevant atomic and molecular features are indicated. 
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